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Models involving large extra spatial dimension(s) have interesting predictions on lepton flavor 
violating processes. We consider some 5D models which are related to neutrino mass generation or 
address the fermion masses hierarchy problem. We study the signatures in low energy experiments 
that can discriminate the different models. The focus is on muon-electron conversion in nuclei, 
/1 —* cy and /j, — > 3e processes and their r counterparts. Their links with the active neutrino mass 
matrix are investigated. We show that in the models we discussed the branching ratio of /1 — > ej 
like rare process is much smaller than the ones of /j. — > 3e like processes. This is in sharp contrast to 
most of the traditional wisdom based on four dimensional gauge models. Moreover, some rare tau 
decays are more promising than the rare muon decays. 

PACS numbers: 13.35.Bv, 13.35.Dx, ll.10.kk 



OO ' !• INTRODUCTION 



In the Standard Model(SM) with fifteen fermions per family neutrinos are strictly massless and the charged leptons' 
weak eigenstates can be chosen to be their mass eigenstates. Thus, each generation has a separately conserved lepton 
number. If one neglects the tiny effects from nonperturbative processes, there is no lepton flavor violating(LFV) 
interaction in SM. However, recent neutrino experiments show strong evidence that neutrinos have none zero masses 
and the three active neutrinos mixp], Q,E|Q. Most physicists take this to be a harbinger of new physics beyond the 
SM. Moreover, finite neutrino masses alone would imply the existence of LFV in charged lepton sector analogous to 
the quarks. If so we expect the Glashow-Iliopoulos-Maiani(GIM) mechanism to be operative in the lepton sector then 
the rate of induced LFV processes will be proportional to the neutrino mass square difference, which is of the order of 
< 10~ 3 (eV) 2 . Hence, they will be hopelessly small for experimental verification. Therefore, additional ingredients are 
' essential for a detectable LFV signature. It is very common in model building to have the new physics that generate 
neutrino masses also given rise to LFV reactions. This link appears to be natural although there is no guarantee that 
O ,' this is case in nature. With this cautionary note we will focus attention to new physics that links the two phenomena. 

Among the numerous beyond the SM models, LFV signatures are most intensely studied in supersymmetric (SUSY) 
ones. The connection with neutrino masses is established through the seesaw mechanism which is the orthodox way 
of getting a small mass for the active neutrinos. Since the latter has a natural setting in grand unified theories (GUT) 
the end result are rather bedecked supersymmetric seesaw models; see e.g. Although the details are different the 
generic source of LFV lies in the mixing of various sfermions. The right-handed Majorana neutrinos play a secondary 



role in this class of models. In general it is natural to expect B(/j, — > ej) ^> _B(/i — > 3e) in SUSY models. 

For non-supersymmetric models neutrino mass generation via the seesaw mechanism would required the right- 
handed neutrinos to be of the GUT scale. In this simplest version all LFV are undetectable. Attempts are now 
made to lower some right-handed neutrinos mandated by the seesaw mechanism to the TeV scale so that the seesaw 
mechanism itself can be tested experimentally. If so then one can optimistically anticipate LFV signatures in the next 
round of experiments [(|. Independent of the details of the models one again expects B(fi — ► ej) B(fi — > 3e) to 
hold true. 

Recently a new avenue has open up in the construction of models beyond the SM that exploits the possible existence 
of extra spatial dimensions. These theories are particularly interesting phenomenologically in the brane world context. 
It is fascinating that many long standing problems in the usual four dimensional (4D) field theories can be overcome 
or take on new perspectives in these higher dimensional constructs. For example the hierarchy problem is solved 
by invoking large extra dimensions. In this note, we would like to draw the readers' attention to the models which 
involve one or more flat extra spatial dimensions. Furthermore, we focus on those that address the neutrino mass 
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problem. In some cases, we predict a reversed pattern of B(/i — > 3e) ^> J3(/i — > e7) compare to SUSY models. On 
the experimental side, it shall be interesting to see this. 

The current experimental limits on muon LFV have already put very stringent constraints on model building. On 
the other hand, the limits from tau LFV arc rather loose. We shall constraint the extra dimension models by the 
muon rare processes data and place upper limits on the rare decays of the r. To avoid any hadronic uncertainties we 
shall focus on purely leptonic processes. We shall also discuss the possible ways to discriminate different models and 
their connections to neutrino masses. 

In this brief review, we give few examples of extra dimension models which give potentially testable LFV signatures. 
These LFV processes are all directly or indirectly related to the generation of neutrino masses. We compare the LFV 
processes in a five dimension(5D) SU(3)w and 5f/(5)@ GUT models where neutrino Majorana masses are generated 
radiatively without a right-handed neutrino which is viable but less discussed alternative to the seesaw mechanism. 
A brief review of this construction is given in |9j . Also, we discuss the LFV processes in split fermion or multi-brane 
scenario. 

The following is our plan for the paper. In sec. II, we will first review the general operator analysis for the lcpton 
flavor violating processes. This will also set the notations for the rest of the discussions. Sec. Ill examines LFV in 
a 5D SU(3) W model. New results of the one-loop calculations are given here. For details of the model and neutrino 
mass generation we refer to pj. In Sec. IV, the LFV processes induced by 5D SU(5) model will be discussed. The 
discussion here has not been presented before. The alternative way of studying the flavor problem using the split 
fermion model is examined in Sec.V. Calculations of the LFV processes in this scenario involves many new unknown 
parameters. The models lack predictive power even semi-quantitatively. However, very general generic trends for 
LFV can be discerned even in this early stages of development. Our conclusions are given in Sec. VI. The necessary 
5D gauge fixing details, which is crucial for loop calculations, are presented in an appendix. 



II. GENERAL OPERATOR ANALYSIS 



First of all, we collect the necessary general formulas for the study of LFV processes. The most important ones are 
the effective interactions of L — I — 7 and L — l — Z where we use the notation L to denote the heavier charged lepton 
which usually is cither /iorr and I is the lighter daughter lepton which can be /i or e. 

In LFV studies, the most important contribution comes from the effective L — I — 7 vertex. The similar vertex 
where a virtual Z replaces the 7 is subdominant in the class of models we are considering. For definitcness we will 
take L = fx and I = e. The most general /1 — e — 7 interaction amplitude allowed by Lorentz and gauge invariance can 
be written as: 

M = -eA;(q)- e (p, - q) j [f EQ (q 2 ) + /mo(<zV] 1» (V - ) + [/mi(<? 2 ) + f E l{q 2 )l b ] u,(p,) (1) 

with the convention e = |e| > used through out this paper and q^ 1 is the photon 4-momentum. For real photon 
emission, only Je\ and /mi contribute. But if a off-shell photon is involved, then all 4 form factors contribute. After 
proper renormalization, the amplitude is finite as q 2 — > 0, so we must have /eo(0) = /mo(0) = 0. It is customary to 
factor out q 2 and rewrite the electric and magnetic form factors as 

f E v(q 2 ) = ^f E o(q 2 ) , fMo(q 2 ) = ^f~Mo(q 2 ) (2) 



and now fEo{q 2 ) and /mo(^ 2 ) finite at q 2 



9 J .C/U \ x l I 1 ,} IV1 U VM J O ■ 



A. L — * hfah and L — > £7 

Using similar notations of the most general effective lagrangian for /i — * 3e and — *■ e'y can be expressed as: 
+ 53(efl7 A 'Mfl)(eR7^e_R) + 54(el7 M Mi)(eL7M e ^) 

+ 55(efl7 M MK)(eT7peL) + 36(eL7 A 'Mi)(efl7M e fl) + h - c - ( 3 ) 
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where 

A « = foitP) + /mi(0)] , A L = -r^j [/mi(0) - f E1 (0)} . (4) 

Also note the anapole form factors f E o and /mo have vector like effective contributions to .93-6 : 

^.93 = Sg 5 = \f E0 (0) - /mo (0)1 (5) 

5.94 = Sg 6 = \f E0 (0) + /mo (0)1 (6) 

which shall be included in the 33,4,5,6- The above effective lagrangian leads to 

B(p^e~f) = 3M7r 2 (\A L \ 2 + \A R \ 2 ) (7) 

flfo->3e) = |gl|2 + |g2|2 +2(| g3 | 2 + | g 4| 2 ) + | g5 | 2 + | g6 | 2 



+64 e a ^ln^-H W|A fl | 2 + |A L | 2 ) (8) 



+8eRe [A R (2g* 4 + g*) + A L (2g* + g* 5 )] 
rrifj, _ 11 
m e 8 

if electron mass is ignored. 

To carry out the calculation, it's convenient to define two dimensionless variables x\ — 2E\jm ll and X2 = 2Eijm lli . 
However, it is important to keep m 2 terms in the intermediate steps in order to properly extract the finite term in 
the last line of Eq.JHJ. Our result agrees with [T(lllll |. 

The expressions of Eq.JSJ), except the last line of dipole operators, can also apply to r — > I + 7 and r — > I1I1I3 
processes . For r — > ee/2, \1\1e processes, the part from dipole operators have double loop suppression from two flavor 
violation vertices and resulting in an insignificant branching ratios thus can be safely ignored. For r decay, the 
branching ratios given above are normalized to B(t — ► ev e v T ). This holds for subsequent discussions of r decays. 

To complete the story, we also give the expression for processes with no identical particles in the final state, namely, 
t — * fj.ee, e/j,fl or h =/= h = h- The above expression for the branching ratio is now modified to: 

B{r^hhh) = |gl|2 + |g2|2 +(| g 3| 2 + |g 4 | 2 + |g5| 2 + |g 6 | 2 ) 
+8ei?e [A R (gl + gl) + A L (g* 3 + g* 5 )} 

+64e 3/ m ^-§W|3 + |A i | 2 ) (9) 
I m 2 2) 

with trivial extension of giS. In arriving the last line of Eq.@, we have ignored the masses difference between m e 
and to m in phase space integration but keep the crucial mass singularity associated with the virtual photon. Not 
surprisingly, the approximation agrees very well with the actual numerical integrations. 

If the photonic dipole operator is the only dominate LFV source, we have the following model independent prediction 
for 

B(t -» e~j) = B{t -> n~j) (10) 
B{r^fiee) _ 2a f ^ to. _ 31 



B(t — > ej) 37r (_ TO e 2 
B(T->efip,) _ 2 «f m ^V _3'. 



B(t — > ery) 37r [ to m 2 
to the accuracy of to 2 /to 2 . To our knowledge, the last two relations have not been presented before 

B. fj, — e conversion in nuclei 

We can write the effective LFV Lagrangian for fj, — e conversion as: 

= e (s - p7 5 ) n 9 (s q - Pql b ) q + ej a (v ~ a7 5 ) \x ^ qj a (v q - a q Y°) q 
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+ \e(t s +t p ^)a a ^^qa aP q + H.c. (13) 

Q 

with self explanatory notations. Here, flavor changing terms in the quark sector are not included since they are not 
expected to be important here. The effective couplings are normalized to {■s/2Gf) 1 ^ 2 ■ For example, the SM Z boson 
has a vector coupling to quarks given by 

v q = T 3 - 2Q sin 2 6 . 

To calculate the conversion rate, we need to promote the interaction from quark level to the nucleon level by 
computing the matrix elements (N\qTq\N) — Gy N NTN with N denotes a nucleon and V = {1, 7 5 , j a) 7q7 5 , <J a p\- 
Since the coherent process is the important one only vector and scalar operators matter: 

(p\q~laq\p) = G q y p pj a p, (n\qj a q\n) = G q y n nj a n (14) 

and 

(p\m\p) = G% V PP , (n\qq\n) = G^nn . (15) 

By conserving of vector current, in the q 2 ~ limit, one can determine that Gy P = G v ,n = 2 and Gy n = Gy P = 1. 
However, one has to rely on the nucleon model to evaluate the scalar operator. For qualitative estimation, we will use 
the result Gs ~ Gy from full non-relativistic quark model but the reader should keep in mind that the uncertainty 
of nucleon model could be as large as few tens percent ^3] • Following the approximations used in [l^, the conversion 
rate, normalized to the normal muon capture rate r capt , can be expressed as |lOt fl3l Il4j : 

£ co „v = PeEeGp ^^ a Z *ff {\4eA L Z +(s- p)S N + (v- a)Q N \ 2 + \4eA R Z + (s+p)S N + (v + a)Q N \ 2 } (16) 

27Hzr capt 

by assuming that the proton and neutron density are equal and the muon wave function does not change very much 
in the nucleus, and F p is a form factor whose definition can be found in 0] and p e (E e ) is the electron momentum 
(energy), E e ~ p e ~ m„. For f 2 Ti( 27 3 Al), F p ~ 0.55(0.66), Z eff ~ 17.61(11.62), and T capture ~ 2.6(0.71) x lO^-^. 
Where the coherent vector and scalar coupling strength of nuclei N are defined as 

S N = s u (2Z + N) + s d (2N + Z), (17) 
Q N = v u (2Z + N) + v d (2N + Z). (18) 

If there are more than one gauge or scalar bosons mediate this process, the above expression can be trivially extended 
with modified couplings: 

M 2 

i Hi 
M 2 

(v±a)Q N ^^2(v i ±a i )Q i NJj f. (20) 

i Zi 

Note that the form factors /eo and /mo in Eq.J2J have extra contribution to the vector couplings: 

2eM w - 2eM w - 2eM w 

= J E0 , da = /mo , ov q = Q q , (21) 

gm p gm p gm^ 

and if Eq.© is the only LFV source, then Ea. (|16fl reduces to the well-known formula given in [l3l | 

8m u F 2 a 5 Zf ff Z , 

B2 onv = eM -{\fMi + Jeo\ 2 + I/mo + Jei\ 2 } ■ (22) 

1 capt 

Also a model-independent relation between the \i — e conversion and the fi — * ej 

ml G 2 F F 2 ^Zj ff Z f \f M i+fEo\ 2 + \fM o + f E i\ 2 
127r3r capt ^ |/mi| 2 + |/bi 



^conv= 10 _^ —( n 12,1/ 12 )S( M ->e + 7 ). (23) 



The above brief review is sufficient for the phenomenological analysis we do. Next, we will head for the extra- 
dimensional models and discuss their LFV signatures. 
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III. 5D SU{3) W UNIFICATION MODEL 

It has been known for a long time that the SM lepton left-handed doublet and the right-handed singlet charged 
lepton in each family can beautifully form an SU (S)w fundamental representation, i.e. L = (e,v,e c ) r [ 0|- This is 
implemented in an electroweak only unification in which SU(2) x U(l) is unified to SU(3)w- One of the attractive 
points of this unification model is the tree level prediction of sin 2 Ow = 1/4. Renormalization group considerations 
point to a relatively low scale of unification at ~ few TeV. We shall use {U ±2 ,V ± } to denote the SU (3)w / {SU (2) x 
U(l)) gauge bosons which have SM quantum number (2, ±3/2). In 4D, the SU(3)w GUT has a fundamental difficulty 
of embedding quarks into SU(3)w representations. This problem can be circumvented by promoting the model into 
five dimensional space time |l7j| and |jj ■ We give a brief summary of the model construction here. 

The extra spatial dimension, with coordinate denoted by y, is compactified into an S\j{Z 2 x Z' 2 ) orbifold. Where 
the circle Si of radius R, or y — [— 7r_R, nR], is orbifolded by a Z2 which identifies points y and —y. The resulting 
space is further divided by a second Z 2 acting on y' = y — ttR/2 to give the final geometry. 

We now have tow parity transformations P : y <-> — y and P' : y' «-> —y' under which the bulk fields can be 
assigned either of the eigenvalues + or -. This freedom is used to break the bulk SU{3)w symmetry to SU{2) x U(l). 
Explicitly, one assigns the following properties to bulk gauge fields 



Mv) = p^(-y)p^ 1 , Mv') = p'A^-y')?'- 1 
My) = -PM-v)?' 1 , Mi/) = -P'AC-yOP'" 1 



(24) 



where the matrices P = diag{+ + +} and P' = diag{+ H — }. Now the (Z%, Z' 2 ) parities of the SM gauge bosons and 
the U, V gauge bosons are (++) and (H — ) respectively. It is easy to work out the Fourier eigenmodes propagating 
in the bulk and see that only fields with (++) parity have zero modes. In other words, only SM gauge bosons have 
zero modes. Both the U, V gauge bosons and all the y— components are heavy KK excitation. Note the second Z' 2 is 
necessary to avoid the presence of zero modes for both SM gauge boson and the exotic U ±2 7 V ± boson at the same 
time. 

The SU(3)w symmetry is explicitly broken to SU(2)l X U(l) at the y = ttR/2 fixed point, where the 4D quarks 
field are forced to live on it. The extra degree of freedom in extra dimensional theories is the key to incorporate SM 
quarks into the SU(3)w symmetry. On the other hand, the lepton fields can be placed anywhere in the bulk or on 
either two fixed points. We choose to put the 4D lepton triplets at y = which is a SU(3)w symmetric fixed point 
so that they enjoy the SU(3)w symmetry. This also avoids possible proton decay contact interactions. 

One Higgs triplet 3 plus one Higgs anti-sextet 6 , denoted as 4>e, with parities 



My) = PM-y) ,My') = P'h(-y') 

My) = PM-y)P-' ,My') = -p' M-y')P'~ l ■ 



(25) 



is the minimal scalar set to give viable charged fermion masses (see |jj )■ Another Higgs triplet 3' with parities (H — ) 
is introduced to transmit lepton number violation essential for generating Majorana neutrino mass through one-loop 
diagrams by a triple Hig gs i nteraction of the type of 3' T 63. This is a 5D realization of radiative neutrino mass 
generation first proposed in [18|, The resulting mass matrix is necessarily of the Majorana type. 
Now we have all the ingredients to write down explicitly the 5D Lagrangian density 



-Tr[G 



MN 



G MN ] 



D-[(D M <p 6 )HD M M] 
+ {D M M\D M M + {D M ^)\D M <t>' z ) 

+ m 




>M 
quark sector 



(26) 



where Gmn, M, N = {0,1,2,3,?/} is the 5D field strength and Dm is the 5D covariant derivative. The cutoff scale 
M* is introduced to make the coupling constants dimensionless. The other notations are self explanatory. The quark 
sector is not relevant now and will be left out. The complicated scalar potential is gauge invariant and orbifold 
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symmetric and will not be specified since it is not needed here. To perform loop calculations, we need to specify the 
5D gauge fixing term, £gf, which will be exhibited later. 

The fields and their parities of this model are summarized below: 

8» = (l,0) ++ + (3,0) ++ + (2,-3/2) + _ + (2,+3/2) + _ 



8 y = (1, 0)__ + (3, 0)— + (2, -3/2)_+ + (2, +3/2)_+ 

bv Ay (u,v)y 

3 = (2,-l/2) ++ + (l,l) + _ 

%i Hs 

3' = (2, -l/2)+_ + (!,!)++ 

fl Wl ff s 

6 = (3,+l) + _ + (2,-l/2) ++ + (l,-2)+_ 

S v ' S v ' S v ' 

where the SM quantum numbers are (SU(2)l, U(1)y) and the subscripts label the parities P, P' . Then it is straightfor- 
ward to obtain the 4D effective interaction by integrating over y and the 4D effective gauge coupling can be identified 
as <?2 = g / \2ttRM* . The orbifold construction is engineered such that there is no tree level LFV in the SM gauge 
interactions. Thus, the success of that model remains intact. But the tree level LFV interaction emerge in the U, V 
gauge interaction which are heavy KK excitation and in the Yukawa interactions. 
The LFV charged current is 



£cc = 92 ^ e L a> i P L (Uiep)ije c Rj U n * + H.c. 

71=1 

+92 £ vn-fPLiUu^ae^V-^ + H.c. (27) 



71=1 



where the subscripts L and R are kept for book keeping. The matrices Ul,r are used to diagonalizc the charged 
lepton mass matrix and Ui ep = U}U R is an extra CKM-like unitary mixing matrix for the lepton sector. 
The LFV Yukawa interactions are given by 

CY = y/2nRM* S Kn [fsiA Vl "j H s,n + f*Hj {^U e L, H% ln + e^Jz^tf^J j) 



71=0 



* 71 = 

V n— 

+ fmj {iW7i e L,i + eL7iZR,j) H W2,n + I ~< R,j :» H il ,) + fsij^RTi e R,j H S2,n + Hx - ( 28 ) 

where K n = (VZ) 1 -^* and 

It = UlfU L , f§ = UlfU R , fp 3 = U T L f^XJ L , /P 3 = Uy^ 3 U L , f H = U R fU L . (29) 
Note that in the new basis the symmetry of fr and fs are not changed. 

A. L — > / + 7 transition 

We begin the discussion by studying a special case that fe^ fa, such that Ur <~ f/£ also /|., and /| are roughly 
diagonal. This hierarchical Yukawa structure is also demanded to yield the observed charged lepton mass hierarchy 
In other words, all the LFV sources are in the Yukawa interaction of <p 3 and (f>' 3 . Since (f>' 3 has nothing to do with 
the charged lepton masses, we can further assume its LFV contribution is larger than <j> 3 , whose coupling is roughly 

~ (m/M w )(f 3 /f 6 ), and f' s 3 ~ 
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In general this class of decays proceeds via the one-loop diagrams. The ones involving the gauge boson U ±2 and 
V ± are suppressed by the GIM mechanism. This leaves the singly charged and neutral scalars as the only possible 
contributors since they both carry two units of lepton charges in the usual scheme. We thus conclude that these 
decays are dominated by the scalar induced Ml and El operators only. Therefore, they provide unique probes of the 
exotic scalar sector. Later we will show that in contrast L — > 31 probes the gauge interactions of the model. 

In this case, the leading contribution loop diagrams are shown in Fig^ Now, briefly discuss the gauge fixing in 



tt'O 
12 Wl 



k \ 



H e (i v, 

(a) 7 n , Z n (b) 

h 



(d) 



(e) 



L e (1 

z ( f ) 



FIG. 1: The leading contributions for n — > e"f and fi — e conversion in the case described in SEC IIII XI The labels n > indicate 
the KK level. 



this model. Because the orbifold parity for <j>' 3 is chosen to be (H — ), it can not develop a VEV and doesn't participate 
the electroweak breaking. The Goldstone bosons consist of the y-components of gauge bosons and the proper linear 
combinations of and <f>§. And the whole 3', H^ 1: H W1 and H' s are physical Higgs. So now it is straightforward to 
carry out loop calculation. For further details, see Appendix. 
The El, Ml form factors are calculated to be: 



f Ll _ 

J mi — 



f Ll 
J El 



384tt 2 MJ 



384tt 2 Mf 



i 



24 



24 



+ If Li) 



(9/w - 7f Li ) 



(30) 
(31) 



where fn — f^ ;i = f^fiL, Mso is the zero mode mass of H<f , and e = (irMsoR) 2 ~ O(0. 1). On arriving at the above 
expression, the contributions of all KK scalar excitation running in the loop have been summed. And if we drop the 
e-terms, the resulting branch ratio can be expressed as: 

2 



B{L^l + ~/) 



967r 3 a 



Ll |2 | I t Ll |2\ 
I J Ml I ) 



El 



768nG 2 F Ml 



E h 



(32) 



= 2.75 x 10" 



/300GeVV 



V M so J 



(fs,le)*fs,eL + (/s,//i)*/s,/iL + (fs,lr)*fs,- 



(33) 



Because the Yukawa couplings of triplet scalars are anti-symmetric, the L — > / + 7 processes have following forms: 



e + j) = 2.75 x 10~ 6 
e + 7) = 2.75 x 10~ 6 



/ 300GeV \ 4 

V Mso J 
/ 300GcV \ 4 

V M so J 



B(r-> Ai + 7) = 2.75 x 10" 



/300GeVY 



V Mso ) 



(fs,eii)*fs,eT 



(34) 
(35) 
(36) 



We have taken M3 = 300GeV as the reference point. If all of the Yukawa couplings are real and none of them vanishes, 
their ratios can be further simplified to: 

B( M - e + 7) : B(r - e + 7) : B(r -> M + 7) = JJ^J ■ JJJ-J ■ jj^j 2 ■ (37) 
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< 2.1x10- 



At this point one can use the data B(/i — > e+7) < 1.2 x 10 11 19] to obtain the constrain (/s 3 eT )*/s 3 /jT 
This is consistent with the expectation from the study of neutrino mass in this as given in 0. There it was found 
that the Yukawa coupling / ^ has to be < 10 -2 and the /'s exhibit the pattern f e ^ > f£. > / ^ Hence it reasonable 
fi — > e + 7 to occur at a rate less than two orders of magnitude below current level. Indeed in the SU(3)w model we 
can link the various L — ► Z7 transition branch ratios to the light neutrino mass matrix elements. Assuming that the 
light neutrino mass are mostly coming from the one-loop quantum correction involving the zero modes of <p 3 and </>6 , 
we have the prediction : 

/ \ 4 

B(fx -> e + 7) : B(t — ► e + 7) : B(t -> fj. + 7) ~ — m 13 m 23 ■ m 12 m 23 : m 12 m 13 (38) 

\ v m T y 

where m,j is the (zj) entry of the light neutrino mass matrix. Interestingly the model naturally accommodates an 
active neutrino mass matrix of the inverted hierarchy type as follows: 

ell 

lee) (39) 
1 e e 2 

where e ~ 0.1. From the above equations, we see that /1 — > ej is suppressed compared to the r — > Z7 decays. This is 
a striking feature of the model. 



B. fi — e conversion 



The fi — e conversion in nuclei will be dominated by the virtual photon exchange. Compared to /i — > e"f it has 
additional contributions from the anapole terms. The corresponding photon E0, MO form factors can be derived as: 



f~Eo(-k 2 ) 



576tt 2 M 2 ^ 

Jyj i=e,p,T 



E 



5767r 2 Af 2 



/ii + + /ii) - ^(/fli + /ii) X/ (2n — 1)2 

G(5 n: Xi) 
^ (2n- l) 2 

n=l v ; 



/ii + — /ii) — ^2 (/fli — /ii) E 



6e 



where <5„ = (—k 2 )/Mf lnl Xi = m 2 /{— k 2 ), i = e,/j,,T, and 



G(<5, a;) = - ln<5 - In a; + - - 4a: + (1 - 2x)yfi + 4xln ^ 4x + 1 — I. 

3 V4^TT+ 1 



(40) 
(41) 

(42) 



As expected, the principal contribution is from the Fig^c) with the zero mode running in the loop. The 
logarithmic enhancements in G(5, x), is due to the exchange of neutral scalars, H$i ( see Fig^a)). Although they 
are suppressed by the KK masses we find them to be compatible to the charged singlet contribution. 
In this process — k 2 ~ to 2 and G has the following limits 



Gen 



m Ji 1 
ln T72 !" o' 



■ In ■ ' l 



1.515, G Tn - - In 



-6.978 



(43) 



The KK sum of these logarithmic enhancements are finite: 



E 



1 



^(2n-l) 

n=l v / H 

So the desired anapole form factors can be expressed as 



tin.. Jl o 

2 In j^f- = — \Ti{m a R) 2 - 0.8362 



(44) 



/bo (m 2 ) = 



/a«)(to 2 ) 



576tt 2 M| 



E 

—e.fj,./ 

E 



5767r 2 Af 2 - 

ou i—e,fi,r 



fu + ^(3/ ffi + /«) + ^|(/ ffi + /L,)[in(TO M i?) 2 + m 

hi - ^(3/fli - /ii) + ^(/fli - /L l )[hi(TO M i?) 2 + m 



(45) 
(46) 
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{VejVuiWr} = { — 1-011,0.837,6.300}. Again, since f s 3 is anti-symmetric, only Jl t — f^ T can contributes. For 
simplicity we assume there is no new CP violation in the scalar sector; then /l = fn and the [i — e conversion rate 
in f§Ti can be expressed as: 



B 1 

conv 



0MB(n 



■7) 



(47) 



if taking 1/i? = 2TeV and Ms = 300GeV as a reference point. It is also possible to have extra contributions from 
KK photon and KK Z excitation, Fig^a-f). One will need to take care of the KK number conservation in the scalar- 
scalar-gauge boson vertices and sum over all the possible combinations. But generally speaking, their contributions 
are further suppressed by (m^R) 2 < 2 x 10 -9 compared to the photon zero mode and can be safely ignored. 

The relation of Ea. H47|) is based on the assumption that Jq 3> fs and <fi' 3 is the dominate LFV source. However, we 
should point out that if f 3 is not so small the neutral scalar zero modes can make /i — > ej and B2 onv compatible and 
deviates a lot from the pure photonic dipole prediction, Eq. I|23|I . 

Again, this demonstrates that L — > Ij and /i — e conversion are very important for us to understand the Yukawa 
structure in the SU(3)w model. 

The question now arises about the photonic dipole and anapole contribution to /i — > 3e. The answer lies in 
Eas. H3QI31l40ll4l|l . We estimated that 



B(fi -> 3e) < 0.04B(> -> ej). 



(48) 



This prediction is not very sensitive to what the Yukawa pattern is. Moreover, the decays L — > 31 have overwhelming 
contribution from other sources of new physics in the model to which we shall turn our attention to next. 

C. L -> 31 

A characteristic of the model is the existence of double charged gauge bosons with LFV couplings. This will induce 
fi — > 3e like processes for the r. In addition there are also KK scalars Ht and Hq which has LFV Yukawa couplings 
which are largely unknown. The Feynman diagrams for the L — > 3/ decays are depicted in Fig|5] Since the Yukawa 



V 



I" 





FIG. 2: Feynman diagrams which lead to r(/i) —* 31 processes. 



coupling are totally unknown, we will postpone the discussion of the contributions from scalars and look at the branch 
ratios, normalized to B(t — > ev T v e ), mediated by U ±2 gauge boson alone first: 

(49) 
(50) 
(51) 
(52) 

(53) 
(54) 



B(t - 


-+3/x) 


= Tx {\U T ^\ 2 




I 2 ) \u»A\ 




B(t 


-3e) 


= Tx (\U Te \ 2 


+ \u eT 


2 )\u ee \\ 




B(t - 


* /zee) 


= Tx (\U T)1 \ 2 




I 2 ) PeeW 




B(t - 


Hfie) 


= Tx (\U Te \ 2 


+ \u eT 


2 ) \u^\\ 




B(t - 


> fiee) 


= |(l^e| 2 + 


\u eT \ 2 ) 


x (|W eAI | 2 4 




B(t - 


enp) 


= + 


l?VI 2 


) x (|Z4 M | 2 - 


f \u^\ 2 



where T = (M W 7ri?) 4 /16 = 1.56 x 10" 5 (2TeV/l/i?) 4 . From the analysis given in sec. II, we know all scalar operators 
give positive contribution. So even though we know nothing about the Yukawa couplings, we can still derive an 
interesting lower bond from the unitarity of Ui ep 



B(r^3e)>Tx \U ee \ 2 (l~\U ee \ 2 ) 



(55) 
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for a given 1/ R. 

If one wants to keep compactification scale 1/R low, say ~ 1.5 TeV, then we would require \U ee \ to be either close 
to zero or one. Furthermore, if we take the upper bound of 1/R < 5 TeV derived from unification seriously we obtain 

B(t -> 3e) > 8.0 x 10- J \U ee \ 2 (l - \U ee \ 2 ) . (56) 

On the other hand, if we assume that the bilepton gauge boson exchange is the dominating FCNC source, another 
interesting upper bond can be derived: 

B(r - 3e)< | = 3.9 x lO" 6 (57) 

with \U ee 1 = 1/ v2 in Eq. Ij55(l . Actually, if all the LFV Yukawa couplings are associated with cf>' 3 as discussed in the 
previous two subsections, the tree-level bilepton scalar contributions to r — > 3e vanish due to the antisymmetry of the 
Yukawa couplings. However, the present experimental limit, 1 — 3 x 10~ 7 [2(j will indicate that the compactification 
radius is closer to the upper limit of 5Tey _1 for this particular case. 



IV. 5D SU{5) MODEL 

The orbifold SU(3)w model discussed above has many interesting and novel features; however, the fact that quarks 
and leptons have to be treated differently is an obstacle towards complete unification. It's a natural attempt to further 
unify the quarks and leptons in a larger GUT group. The simplest group for that is SU(5). Now all fcrmions are on 
equal footing and can be clustered into 2 SU(5) representations, i.e. \&g = {d c ,L},^ 10 — {Q,u c ,e c }. 

Similar to the SU(3)w model, the model is embedded in the background geometry of Si/Z 2 x Z 2 orbifold. The 
bulk SU(5) gauge symmetry is broken to the SM by orbifold parities , with parity matrices diag{+ + + + +} and 
diag{ h +} for Z 2 and Z 2 transformations respectively. These are generalizations of the SU(3)w case. 

Since no right-handed neutrinos are added, neutrino masses can be generated through quantum correction by using 
either 10 or 15 bulk scalars plus the 5'(10/15)5 interaction mandated by breaking to the SM gauge group. The 
orbifold parities of 10 or 15 bulk scalars are determined to be (++) by considerations of proton decay. They split 
into following components: 

15 S (++)=P 15 (m,-~J +T 15 (l,3,l) ++ + Ci 5 (s,2,~J , 



10„ 



Pw 3,1, 



Sio(l,l,l)++ + Ci 3,2, 



A careful analysis shows that by using 15(10) the resultant neutrino mass matrix favor the normal(inverted) hierarchy 
8J. It was also found that extra fine tuning efforts were needed to obtain phenomenologically acceptable neutrino 
mass patterns by using 10 alone; so we will only discuss the case which implements 15. 

The P components induce tree-level K° — K° mixing. To satisfy the experimental constraints, it is required that 
Mp > 10 5 GeV. On the other hand, the two bulk Higgs in 5, 5' which are responsible for the SM clcctrowcak symmetry 
breaking share the same (-1 — h) parities as 15. The brane Yukawa interaction term is easily constructed to be 



15 



{AB} 
15 



H.c. 



(58) 



where A, B are the SU (5) symmetry indices. It can be seen that to contain the necessary LFV source to generate 
neutrino Majorana masses. The neutrino mass matrix elements are proportional to {M.)^ oc m fc/ifc//fc where k 
are the generation indices and uik is the mass of /c-charged lepton running in the loop. 

The extra Higgs doublet in the 5' is good for gauge unification. By adding additional decaplet bulk fermion pair 
with (-1 — ) parity and mass around 10 — 120 TeV, the unification is achieved at 3 x 10 16 — 10 15 GeV or equivalcntly 
1/R ~ 10 14 GeV. The high scale unification or tiny radius of extra dimension makes KK excitation decouple from 
most phenomenological studies and basically we only need to consider the zero modes. 

Below unification scale or equivalcntly the low energy 4D effective theory is a two Higgs doublets like model. In 
general the two Yukawa patterns are not aligned which that can lead to severe tree level charged neutral flavor 
changing (FCNC) interaction. A Z 2 symmetry is usually assumed to forbid such tree level FCNC (2L| - In this model, 
there is no such freedom since the Yukawa patterns are determined by the geometrical setup. The ^iq of the first 
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two generations are assigned to be bulk fields and the other fermion fields, and ^g' 2 ' 3 , are localized at the y = 
brane. In doing so, the salient SU(5) prediction of mb/m T ratio is preserved and give small hierarchy patterns in the 
Yukawa couplings of both 5, 5' scalars, i.e. 

/ 8 8 1 \ / S 2 5 2 5 \ 

Vdoc \ 8 8 1 ,!/»« \6 2 S 2 8 . (59) 
\8 8 1 J \ 8 8 1 J 

Due to volume dilution factor we get 8 ~ 0.1 which measures the amount of overlap between brane and bulk fields. 
The specific Yukawa pattern above successfully generates mass and mixing hierarchy of charged fermions: 

mi, : m s : rrid — m T : : m e ~ 1 : 8 : 8 2 (60) 

m t : m c : m u ~ 1 : 8 2 : S 4 , (V us , V cb , V ub ) ~ ((5, <5, 5 2 ) (61) 

The rotation from weak to mass eigenbasis simultaneously diagonalizes the two Higgs doublets Yukawa couplings. 
Thus, we do not have the FCNC problem due to mixing between two Higgs doublets. 

Instead, now tree level LFV processes can be mediated by the triplet component T15 in 15. The only important 
ones are the \i — > 3e like processes, see Fig|3] 



d s 



P15 



Pv 



d d 



s d 



1 I 



FIG. 3: Feynman diagrams for tree level K — K mixing and /i — ► 3e process. 



An explicit calculation give the branching ratio of fi — > 3e: 



^^.ffltf^l . (62) 



g%(nRM*) 2 \ M T 

The mass difference AM^ in K° — K° mixing arises from P10.15 can be used to eliminate the ambiguity of absolute 
strength of Yukawa couplings. The ratio of Yukawa couplings can be replaced by the ratio of the corresponding 
elements in M. v . Since only the 15 Higgs is used , we have 

16 / AmgA 2 / M P \ 4 ( 2m 11 m 12 



Br(n -> 3e) ~ 3.02 x 10- ib — — - x ±i££ . (63) 

V ; \Am K J \Mt) \m lim22 + (2^m 12 ) 2 J y ' 

It is straightforward to extend the analysis to r — * 3/ transitions. Assuming that the hierarchy of the elements of 
neutrino mass matrix is smaller than factor 100, this model predicts 



Br{n — > 3e) : Br(r — > 3e) : Br{r — > 3/i) : _Br(r — > \iee) : Br(r — > e/x/x) 

\ 4 2 / \ 4 o / \ 4 o / \4 9 

m M \ mj 3 / m e \ m 23 / to m \ m 23 / m e \ m\ 2 



TO 22 \ m T J ttl 22 \m T J mi! \m T / V 77 ^/ m\\m 22 

The photonic form factors due to T ±2 , T scalar one loop diagrams can be obtained: 



(64) 



2 c/-15/fl5\* 

" 16^ 2 Af 2 24 ^ 

2 /15/fl5\* 
JmO-JeQ ~ 167r 2 M 2 g 



G(<5, Xj) - i 



(66) 



The chiral structure in the above result is easily understood because that only the lepton doublets interact with triplet 
scalar. 

Due to the factor (m^/My) 4 , the fi — > ey process is strongly suppressed. Taking Mt — 10 5 GeV, the \iTi — > eTi 
conversion rate is estimated to be ~ 1.2 x 10 1 4 1 y (// e 5 )* | 2 where we have kept only logarithmic terms which is 
sufficient for an order of magnitude estimate. The experimental bound is 3.6 x 10~ n 0] which is not stringent. On 
the other hand the recently proposed experiment aimed at a detecting a signal at the 10 -16 — 10 -17 will be very 
encouraging. Even a negative result will provide stringent constrain on the otherwise unknown Yukawa couplings. 
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V. SPLIT FERMION MODEL 



An interesting scenario was introduced by |22( to solve the charged fermion masses hierarchy problem. The basic 
idea is to postulate that fermions are bulk fields and they interact with a non-dynamic background scalar potential. 
In the 5D version the bulk fermions are vectorlike, but only one of the chiral zero modes will be localized at the zero 
of the background potential modulated by the 5D mass terms. The chirality of the zero mode is determined by the 
sign of slope of the background potential at the zero point. The fermion zero modes are given a Gaussian profile in 
the fifth dimension and each has its own unique position in the extra dimension. The widths are controlled by the 
potential slope at the localized position. For simplicity, we will assume a universal width for all the SM fermions. The 
5D fermion excitation are vectorlike and will be located at the same position of their zero modes. Roughly speaking, 
the energy gap is ~ l/(Width) 3> 1/R 3> Mw and they decouple from the phenomenology we are interested in. Since 
the SM fermions are scattered over the fifth dimension, to preserve the gauge interaction universality the SM gauge 
fields are forced to be bulk fields too. To illustrate the basic physics involved it suffices to build a model on an S1/Z2 
orbifold so that one can remove the unwanted y-components of gauge bosons zero modes which are identified with 
the SM gauge bosons. However, to break electroweak symmetry, a dynamic bulk Higgs is necessary. 

To be more concrete we take the extra dimension to be the interval y G [— irR, -kR] and the fermions are fixed in 
different positions Zi in this interval. The usual Kaluza-Klein ansatz is invoked that any bulk field factorizes into a 
4D field times a 5D wave function. Thus, for the fermion field we have ipi(x,y) = g(zi,y)ip(x) where g(zi,y) is taken 
to be Gaussian distribution in the fifth dimension: 



l(zi,y) 



1 



(7To£)V4 



exp 



(y - *i? 

2„% 



where <jq is the universal width of Gaussian distribution. If <jq <C R, g acts like the Dirac delta function, effective 4D 
interactions are obtained by integrating out the y direction. Then any pair fermions get an exponential suppression 



g( z i,y)g( z 2,y) = exp 



(zi - z 2 f 



z\ + z 2 
9 I o 'V 



as a result of integrating over Gaussian functions. Therefore, the linear displacement between left-handed and right- 
handed fermions in the fifth dimension translates into exponential Yukawa hierarchy in 4D theory. One set of solution 
for the quarks positions have been found numerically that can accommodate the mass hierarchy and the CKM 
mixing 23] . To accommodate the CP violation, the overall Yukawa coupling strength of up and down type quarks 
must be different |24l l25| . Although a fundamental theory of where to place the fermions are located is lacking, we 
have at least one realistic solution for where the quarks are located in the extra dimension. For the lepton sector, we 
don't have enough constraints to pin down the solution. As pointed out bv pHI^ . the FCNC interaction is induced 
geometrically where phenomenological constraint in quark sector have also been discussed. In fact, flavor violation is 
generic in any multi-position models. This comes from the fact that weak to mass eigenbasis rotations can only make 
the SM interactions (zero modes) flavor diagonal, the KK modes cannot be simultaneously diagonalized. 

In general, the LFV can be discussed independently of the neutrino sector. The 4D effective LFV lagrangian can 
be expressed as 



C 



LFV 



E F!) *f 



9hUn{ij)L + g R U n ^R IjZ™ 



U n(ij) L + U n(ij) R 



E 



V2 



kru^L^w;- + h.c. 



(67) 



where A n , Z n and W n are the n-th KK excitation of the photon, Z and W bosons and gt,lR — ^MO — Qi sin 2 9w- 



L I R 

The matrices U n are a combination of the unitary transformations Vl ir that take the lepton weak eigenstates to 
their mass eigenstates and the cosine weighting of the n-th KK modes. Explicitly, they are 



U n /R = ^diag ( cos 



L/R 



R 



L/R L/R" 

ny 9 nyn 
cos — - — , cos ■ 



R 



R 



V, 



L/R 



(68) 



for n = all the cosine factors become one and the interactions reduce to SM as expected. It is clear that the JJ L I R 
is no more diagonal nor unitary in general. 
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1 

Q 

FIG. 4: Flavor changing due to KK gauge bosons. 




In this model, the r — ► 3Z decay and fj, — e conversion happen at tree-level, see Fig^J The corresponding effective 
LFV couplings are: 



.93 = (2A%i?) 2 (sin 2 8 W + g%/ cos 2 6 W ) £ U ^<^ ( 69 ) 

71 = 1 

.94 = (2A%fl) 2 (sin 2 9 W + g|/ cos 2 fl w ) V ^ffi" (70) 

n=l 

.95 = (2M w R) 2 (sm 2 W + g fl g L / cos 2 9 W ) V " ,et 2 (71) 

n=l 

.96 = (2M w i?) 2 (sm 2 g w + g L9R / cos 2 M V -S-^jg (72) 

n=l 

and the a, v q and a q can be obtained in a similar way. Also the lepton universality is broken due to flavor dependent 
couplings in the KK gauge interaction. We refer the reader to |26j for a detailed analysis. 

The leading /i — > ej contribution comes from the one loop corrections. We need to fix the gauge before proceeding. 
The necessary details of 5D gauge fixing are collected in the appendix. After properly identifying the Goldstone 
boson, the usual 4D gauge technique can be straightforwardly applied here. Note that the Yukawa couplings of 
the physical KK scalars are suppressed by the factor of ~ (miR)/n. Although there is residual GIM cancellation in 
the KK gauge boson interaction, we expect the leading LFV are from KK gauge interaction. 

The LFV photonic form factors due to KK gauge boson and their Goldstone boson can be calculated. The KK W 
bosons' contribution to the photonic form factors are given by: 



,11 _ , il _ 7 gg \ " i m ^R) 2 TT L* T jL f 7 o\ 
/Ml - J El - ^[75^2 ^2 U n,ie u n,iiJ,> 

~ n—1 



fW _ fW _ 23 92 \ - ( m ti R ) 2 t T L* t T L ( 7A \ 
JMO - -lEQ - 72 iQ n 2 n 2 U n,ie u n,in K'^) 

n—1 



and for the KK Z bosons they are 

\2 r ,2 



R - 1 



z _ (m n R) 2 g\ \ 1 f 1 \2jjL* jjL + ri 2jjR*jj 

J m\ I ei - 87r 2 cos 2 6 ^ ^ n? \ 3^% e >v ^«>"=^ 

n—1 i 

+ — L 9LgR [Un* ie Un,in ± U n*ie U n,i^ \ > 
"V ) 

r EO/ Mo = Jm 2 $\oL £ £ h { ( G{Sn > Xi) + 5) ^ieKi, ± TO } ( 76 ) 

" ' n—1 i ^ ^ ~ ■ 



(75) 



Similar contributions from KK photons can be easily read from the above by replacing (52/ cos0) — * e, gi, — > 1, and 
9R -> 1- 

These photonic form factors give extra contribution to fi — > 3e and y, — e conversion processes but can't compete 
with those tree-level KK gauge boson exchanging diagrams. However they are the sole sources of new physics for the 
L — > I + 7 process. 

In addition to the usual ignorance with regard to Yukawa coupling there are more unknowns in the lepton locations 
and the Gaussian widths. Ad hoc simplifying assumptions have to be made. Hence, this kind of model suffers from a 
lack of predictive power in LFV studies. More data such as the scale of neutrino mass and more complete knowledge 
of the neutrino mixing matrix will help greatly. However, we can extract some generic features for this kind of models 
as follow: 



13 



1. fi/r — > 3? and /xTi — > eTi(or t — > Z + hadrons) will happen at tree-level from the exchange of KK scalars, 
photons and Z bosons. 

2. L — > Z7 proceeds at the one-loop level and hence is expected to be suppressed compared to the previous modes 

3. Violation of lepton universality will occur. The best signal will be to look for the violation in W — > Zji/j decays 

26]. 

Unfortunately a more quantitative statement about the level of the effect eludes us for now. 

VI. CONCLUSION 

We have studied and reviewed LFV processes in 5D gauge models that are related to neutrino mass generation or 
address the flavor problem. Specifically we focus on two complete models which generate neutrino masses radiatively. 
This allows us to see in detail how the two issues can be related. The models are based on SU(3)w and SU(5) 5D 
unification. They give rise to different neutrino mass patterns 0,13; thus, it is not surprising that they give different 
prediction for LFV. The SU(3)w model has a unification scale at ~ TeU and makes essential use of bileptonic scalars. 
It also contains characteristic doubly charged gauge bosons. The SU (5) model is a 5D orbifold version of the usual 
GUT. The unification scale is much higher at 10 15 GeV. The important ingredient for LFV and neutrino masses is 
the 15 Higgs representation. The triplet Higgs of this model plays the crucial role here. 

We found that for the SU(3)w model the rare r decays are much more enhanced compare to their counterpart /j, 
decays. Among the r — > I + 7 decays the largest mode is the [i + 7. Even for this mode we expect it to be < 10 -14 
which is much lower then current experimental reach. 

The decay modes r — > 31 have a better chance of being observe. This stems from the fact that they are tree level 
processes induced by the bileptonic gauge bosons or scalars. Since they are KK modes they have high masses controlled 
by the extra dimension compactification radius which is < 5TeU from consistency and unification considerations. An 
order of magnitude improvement on the current limit will be valuable information on the unknown Yukawa couplings. 

For the orbifold 5D SU (5) model the muon to electron conversion in nuclei can be within the experimental capability 
of the proposed experiment at Brookhaven National Laboratory!^. As in the previous model fi — > e + 7 will not be 
observable. This is very different conventional 4D unificational models. 

The split fermion model also have the characteristic of L — * 3/ and fj, — ► e conversion dominating over L — ► Z7. 
We cannot be more quantitative due to proliferation of unknown parameters. This model have lepton universality 
violation which is not present in the previous two models. This can serve as a differentiating tool. 

It is clear that in order to unravel the physics behind the flavor problem all modes of LFV must be searched for. The 
usual 4D supersymmetric model will favor L — > Z7 where as the 5D models prefer L — > 31 and/or /1 — * e conversion. 
To this we add lepton universality test as a probe. 
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APPENDIX A: 577(2) x U(l) IN 5D S1/Z2 MODEL WITH A BRANE AT y = 

Now we present the gauge fixing scheme used for the 5D electroweak interaction with one bulk Higgs doublet. For 
simplicity the background geometry is S\/Z2- The fifth gamma matrix was chosen to be 7 a = 17 s . The 5D Lagrangian 
is 

C 5 = -1 F MN F MN - \g^ mn G m ] n + p„$) f (D M ^) +■■■ (Al) 

where 

Fmn = 9m Bn 

MN — °MA N 



- d N B 



n-°m ■ 



d N A[ 



92 abc Ab) Ac) 
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B and A stand for the U(l) hyper charge and SU(2) gauge fields respectively. In this convention, Q = T3 + Y. We 
adopt the usual conventions: = ^(^f T iAf), P M (hoton) = (c w B M + s w Af) and Z M = {c w Af - s w B M ), 

or B M = c w P M -s w Z M , Af = c w Z M + s w P M ), where c w = h/VW+M and s w = gx/VW+H- 9s = VW+M 
are introduced to simplify the notation. The symmetry breaking pattern is same as in the usual 4D SM. The bulk 
Higgs doublet acquires a nonzero VEV after SSB, 

* = f <h+h° ) ,h Q = Q + iX °. (A2) 
The generalized linear gauge fixing is introduced (28] (other schemes can be found in 29]), 



2 



- j\d^W+^ + t[d y W+> v -iM w h+]\ 2 

- ±(d,z» + z[d y zy-M zx a ]) 2 . (A3) 

to remove the mixing between gauge bosons and Higgs. Therefore, the Goldstone bosons and the physical scalars can 
be easily identified: 

Gl = py , AfW = n /R (A4) 

G n [ c ri s n Xn] ' \_ s n c n Xn\ 

MP = ^Jn 2 /R 2 + M 2 , 4 = Mz/M™ (A5) 

M&> = ^/n 2 /R 2 + M 2 v , s z n = M W /M$ (A6) 
H° = 00 , 4'' = ^JrF+MI (A7) 

where G° and G ± are the KK Goldstone bosons, S° is the physical KK pseudo scalar, and H°,H ± are the physical 
KK scalars. The usual R$ gauge can be extended to the 5D S1/Z2 model with little modification, like My/ M^ 
and so on. 

The Goldstone bosons are mainly constituted by the fifth gauge components with a small fraction of KK Higgs 
bosons mixed interaction. On the other hand the Goldstone bosons couple to brane fermions through their Higgs 
components. In contrast the physical scalars are mainly composed of KK Higgs plus small amount of the fifth 
components of gauge fields. 

This scheme can also be applied to the models built on the S\j{Zi X Z' 2 ) orbifold with little modification. 



APPENDIX B: GAUGE FIXING FOR THE ORBIFOLD MODELS WITH MORE THAN ONE SCALARS. 

The method can be easily extended to the cases with multi scalars. Taking a 5D two Higgs doublets Model(2HDM) 
as an example, with VEVs (tpi) = v\, (<fo) = V2 and tan/3 = V2/V1, the physical charged scalars and pseudoscalars 
are 

H ± = sin 0<f>f - cos (3$ , A = sin filmed - cos /3im0°, (Bl) 

just like the usual 4D 2HDM. The only difference is that the orthogonal linear combinations a — cos 0Im4>\ + 
sin /3Jm^2 and g ± = cos /3(j)f + sin /3<fif will mix with the fifth components of gauge fields to form the real Goldstone 
bosons: 

G° = cos 9°V° " sin 6»°a° , sin 0° n = Mq/\J M§ + n 2 /R 2 , (B2) 
Gt = cos 6i V± T i sin 6±g± , sin fl± = My/ \J M v + n 2 / R 2 . (B3) 
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